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ABSTRACT 

We present multiwavelength observations for a large sample of microjansky radio sources detected 
in ultradeep 1.4 GHz maps centered on the Hubble Deep Field-North (HDF-N) and the Hawaii Survey 
Fields SSA13 and SSA22. Our spectroscopic redshifts for 169 radio sources reveal a flat median redshift 
distribution, and these sources are hosted by similarly luminous optical L* galaxies, regardless of redshift. 

We suggest that the absence of low optical luminosity galaxies at low redshifts, where there are no 
selection effects, is due to small galaxies not being efficient at retaining the cosmic rays necessary to host 
microjansky radio sources. This is a serious concern for radio estimates of the local star formation rate 
density, as a substantial fraction of the ultraviolet luminosity density is generated by sub-L* galaxies at 
low redshifts. From our submillimeter measurements for 278 radio sources, we find error-weighted mean 
850p,m fluxes of 1.72±0.09mJy for the total sample, 2.37±0.13mJy for the optically-faint (/ > 23.5) 
subsample, and 1.04±0.13mJy for the optically-bright (/ < 23.5) subsample. We significantly (> 3cr) 
detect in the submillimeter 50 of the radio sources, 38 with / > 23.5. Spectroscopic redshifts for three 
of the I < 23.5 submillimeter-detected radio sources are in the range z = 1.0-3.4, and all show AGN 
signatures. Using only the submillimeter mapped regions we find that 69 ± 9% of the submillimeter- 
detected radio population are at / > 23.5. We also find that 66±7% of the Ss 5 o Aim> 5 mJy (> dtr) sources 
are radio-identified. We use our spectroscopic sample to determine the evolution with redshift of the 
radio power, and hence the far-infrared (FIR) luminosity (through the local FIR-radio correlation). We 
find that millimetric redshift estimates at low redshifts are best made with a FIR template intermediate 
between a Milky Way type galaxy and a starburst galaxy, and at high redshifts with an Arp 220 template. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: formation — galaxies: 
starburst 

1. INTRODUCTION 

Some of the deepest 1.4 GHz maps ever obtained with 
the Very Large Array (VLA) are centered on the Hub¬ 
ble Deep Field-North (HDF-N) and the Hawaii Survey 
Fields^^ SSA13 and SSA22. These maps should be sen¬ 
sitive to the synchrotron emitting disks of galaxies out to 
z ~ 1 and the nuclear starbursts of ultraluminous infrared 
galaxies (Sanders & Mirabel 1996) out to z ~ 3. Most 
importantly, since 1.4 GHz emission does not suffer from 
absorption by galactic dust or neutral hydrogen at high 
redshifts, even the most distant, dust-obscured galaxies 
should be transparent. The 1.4 GHz maps of the HDF- 
N, SSA13, and SSA22 fields are not only unprecedented 
in their sensitivities (e.g., the r.m.s. of the SSA13 image 
is 5 /rJy) but also in their detailed resolution, which has 
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enabled unambiguous optical identifications to be made 
on images aligned to within 0.1" — 0.2" of the radio FK5 
coordinate grid. These optical identifications are com¬ 
posed of 10%—20% late-type Seyfert galaxies and bright 
field elliptical galaxies and 60%—70% moderate redshift 
(0.1 < z < 1.3) starburst galaxies characterized by their 
diffuse synchrotron emission on galactic scales (1" — 2"; 
T. W. B. Muxlow et al., in preparation). The remain¬ 
der of the sources are unidentified to I ~ 25 Thus, the 
faint radio source population appears to be dominated by 
moderate-redshift star-forming galaxies (e.g., Windhorst 
et al. 1995; Richards 2000). 

Sensitive submillimeter follow-up observations with the 
Submillimetre Gommon User Bolometer Array (SGUBA; 
Holland et al. 1999) on the 15 m James Glerk Maxwell 
Telescope^^ of the optically-faint tail of the radio source 
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population have revealed significant numbers of 5'850/im> 
5mJy sources (Barger, Cowie, & Richards 2000, hereafter 
BCROO; Chapman et al. 2001, hereafter COl; Chapman 
et al. 2002a, hereafter C02) whose submillimeter emission 
appears to be produced primarily by star formation rather 
than by accretion onto supermassive black holes (Barger 
et al. 2001a, b). 

A radio selection technique can therefore be used to 
identify bright submillimeter sources without the need 
for long submillimeter exposures. COl found that high- 
redshift S'850;im>5mJy sources can be detected with 
SCUBA in photometry mode at a rate of about one per 
hour, which is significantly more rapid than the rate of 
about one per 10 hours achieved in blank-field jiggle map 
surveys (e.g., Barger et al. 1998; Hughes et al. 1998). 
Moreover, the S' 850 /im/<S'i. 4 GHz ratio can be used as a mil- 
limetric redshift indicator (Carilli & Yun 1999, 2000, here¬ 
after CY99 and CYOO; BCROO; Dunne, Clements, & Bales 
2000, hereafter DCEOO). Millimetric redshift estimates of 
submillimeter-detected radio sources are typically in the 
range 1 < z < 3 (e.g., BCROO; Smail et al. 2000; COl; 
C02; Ivison et al. 2002). 

Flux-limited radio selection does introduce a luminos¬ 
ity bias with redshift (see § 5.1); even very sensitive radio 
maps, such as those used in the present paper, will not de¬ 
tect galaxies at very high redshifts. C02 modeled the radio 
selection function in detail using different galaxy evolution 
scenarios and found that the primary bias was a cut-off of 
galaxies at z < 3 due to the radio sensitivity, although the 
specific form of galaxy evolution and dust properties also 
affected the radio-detectable sources. However, the large 
percentage (~ 70%; see § 4.3) of S'85o^m>5mJy blank- 
field counts recovered through radio selection shows that 
most of the bright submillimeter sources can be found this 
way, in agreement with predictions from the C02 best-fit 
models. 

Although micro) ansky radio sources have proved to be 
a crucial resource for identifying and studying the nature 
of submillimeter-luminous galaxies, the submillimeter- 
detected samples are only a small fraction of the radio 
population. In this paper we present multiwavelength ob¬ 
servations of a large sample of microjansky radio sources 
detected in ultradeep 1.4 GHz maps of the HDF-N, SSA13, 
and SSA22 fields to learn about the nature and evolution 
with redshift of the microjansky radio population, as well 
as the implications of using this population to study the 
bright submillimeter source population. We assume a flat 
Ha = 2/3 Universe with Hq = 65 km s“^ Mpc“^. 

2. OBSERVATIONS 
2.1. Radio imaging 

Data at 1.4 GHz were obtained for the HDF-N, SSA13, 
and SSA22 fields using the VLA. The HDF-N radio map 
(Richards 2000) was observed in the A-array configuration, 
resulting in an ~ 1.5" beam and reaching a sensitivity of 
40 ^Jy (5 (t). The SSA13 radio map (E. B. Fomalont et 
al., in preparation) was observed in both the A and B- 
array configurations for a total of 100 hours, resulting in 


an ~ 1.5" beam and reaching a sensitivity of 25 /iJy (5 (t). 
The SSA22 radio map (M. Yun et al., in preparation) was 
observed in the B-array configuration for 12 hours, result¬ 
ing in an ~ 5" beam and reaching a sensitivity of 60 /iJy 
(5cr). 

2.2. Optical and near-infrared imaging 

Deep /-band imaging was used to identify the radio 
sources. The /-band image of the HDF-N field was ob¬ 
tained with the Subaru Prime Focus Camera (Suprime- 
Cam; Miyazaki et al. 1998) on the 8.2 m Subaru^^ tele¬ 
scope. Details of the observations and reductions can be 
found in Barger et al. (2002), and catalogs of all the galax¬ 
ies and stars in the field are presented in P. Capak et al., 
in preparation. The HDF-N optical image covers 24' x 24' 
to / = 25.8 (5 (t). In the HDF-N field, 357 radio sources 
have /-band measurements. 

The /-band image of the SSA13 field (Barger et al. 
2001a) was formed from a mosaic of overlapping images 
obtained with the Low Resolution Imaging Spectrograph 
(LRIS; Oke et al. 1995) on the 10 m Keck^"* telecopes. In 
the SSA13 field, 145 radio sources have /-band measure¬ 
ments. 

The /-band image of SSA22 was obtained with the 
CFH12K camera (Cuillandre et al. 2000) on the 3.6 m 
Canada-France-Hawaii Telescope (CFHT). The 8.1 hour 
integration reaches I = 26.1 (5(t). Details are presented in 
S. C. Chapman et al., in preparation. In the SSA22 field, 
122 radio sources have /-band measurements. 

For / > 21 we used 3" diameter aperture measurements 
corrected to total magnitudes using an average median off¬ 
set. For / < 21 we used isophotal magnitudes measured to 
1% of the peak surface brightness of each galaxy. Isophotal 
magnitudes provide a better measure of the total magni¬ 
tudes in these generally low redshift galaxies with extended 
angular sizes. 

Deep //-band imaging of SSA22 was obtained at UKIRT 
using the UKIRT Fast Track Imager (UFTI). The small 
UFTI field (50" x 50") was centered on the positions of the 
> 3(7 submillimeter-detected radio sources. Each source 
was imaged for a total of 3600 s using individual exposures 
of 60 s; the limiting magnitude in a 2" diameter aperture 
is K = 20.8 (5(t). The fast tip/tilt, adaptively-corrected 
imaging resulted in seeing better than the median condi¬ 
tions of 0.45" FWHM. The data were reduced using the 
Starlink UKIRT/UFTI image processing tools under the 
ORACDR environment (Bridger et al. 2000). We wrote cus¬ 
tom ORACDR scripts to Optimize point-source sensitivity in 
our essentially blank-field observations, creating flat-fields 
from each 9-point dither and high signal-to-noise thermal 
background images from 60 minutes of data. 

For sources in the HDF-N, //-band magnitudes were in¬ 
ferred from the deep HK' image of the Hubble Flanking 
Fields (Barger et al. 1999; HK' — K = 0.3) obtained using 
the University of Hawaii Quick Infrared Camera (QUIRC; 
Hodapp et al. 1996) on the 2.2 m University of Hawaii 
telescope and the CFHT. The image reaches an effective 
limiting magnitude of // = 20.1 (5 (t). 
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Fig. 1.— Configuration of the submillimeter data in the HDF-N, SSA13, and SSA22 fields (as indicated in panels a,b,c). A central area of 
~ 120 arcmin^ has been mapped with SCUBA to an uniform depth of ~ 2 mjy in all three fields, although small subregions go much deeper. 
For illustration, we have indicated the radio sources with squares within the field with the deepest radio map (SSA13). Targeted submillimeter 
follow-up of 7 > 23.5 radio sources was undertaken outside of the SCUBA mapped regions, but within the VLA primary beam (shown with 
filled circles), as well as within the SCUBA mapped regions to confirm map detections. The half-power point of the VLA primary beam at 
1.4 GHz is 30' diameter. 


2.3. Optical spectroscopy 

Spectroscopic measurements of 158 radio sources in the 
HDF-N and 51 in the SSA13 held were made with LRIS 
on Keck and with the HYDRA spectrograph (Barden et 
al. 1994) on the 3.5 m WIYN telescope. We hereafter re¬ 
fer to the radio sources in the HDF-N and SSA13 fields 
with spectroscopic measurements as our spectroscopically- 
observed radio sample. Redshift identihcations were ob¬ 
tained for 169 sources: 126 in the HDF-N and 43 in SSA13. 
The spectral observations and reductions are presented 
and described in BCROO and E. B. Fomalont et ah, in 
preparation (see also Cohen et al. 2000 and Barger et al. 
2002 for many of the redshifts in the HDF-N held). 

2.4. Submillimeter imaging 

The submillimeter data were taken in SCUBA jiggle 
map, SCUBA raster map, or SCUBA photometry mode 
at 850 /imduring a number of observing runs with good to 
excellent conditions. The conhguration of the data in all 
three helds is similar: a roughly 10' x 10' region in the cen¬ 
ter of the 40' diameter VLA primary beam was mapped 
with SCUBA in either jiggle map or raster map mode to an 
average sensitivity of 2 mJy r.m.s. (with small subregions 
reaching as deep as 0.5mJy r.m.s.). These maps present 
an unbiased study of the submillimeter properties of the 
microjansky radio sources. The SCUBA jiggle maps are 
described in detail in Barger, Cowie, & Sanders (1999), 
BCROO, and Barger et al. (2001a, b), and the SCUBA 
raster maps are described in Borys et al. (2002) and C. 
Borys et al., in preparation. For the jiggle and raster 
maps, submillimeter huxes were measured at the posi¬ 
tions of the 1.4 GHz (> 5cr) sources from the catalogs of 
Richards (2000; HDF-N), E. B. Fomalont et al., in prepa¬ 
ration (SSA13), and S. C. Chapman et al., in preparation 
(SSA22) using beam-weighted extraction routines that in¬ 
clude both the positive and negative portions of the beam 
profile. 

Additional targeted SCUBA photometry mode obser¬ 
vations were made of the majority of the optically-faint 


radio sources over a wider part of the VLA primary 
beam, including some sources in regions already mapped 
by SCUBA. The photometry data are described in COl, 
C02, and S. C. Chapman et al., in preparation. These 
observations were centered on the radio source coordi¬ 
nates producing a flux measurement at that position. 
The configuration of the submillimeter data in each of 
the three fields is depicted in Fig. 1. Of the 278 radio 
sources in the three fields with submillimeter measure¬ 
ments, all have 850 /rmflux measurements with uncertain¬ 
ties less than 5mJy, and of these, 219 have uncertainties 
less than 2.5 mJy. In total, we signihcantly (> 3cr) detect 
in the submillimeter 50 of the 278 radio sources. Twenty- 
eight of these lie in the jiggle or raster maps, and 22 
are from photometry mode. For reference, there is about 
one blank sky SCUBA source with > 5 mJy for every 200 
SCUBA beams. 


3. OPTICAL PROPERTIES AND SPECTROSCOPIC 
REDSHIFT DISTRIBUTION OF THE MICROJANSKY 
RADIO SOURCES 

Richards et al. (1999) presented the overall I magnitude 
distribution for the 1.4 GHz and 8.4 GHz samples in the 
HDF-N and the 8.4 GHz sample in the SSA13 held. They 
found that the optical counterparts had a wide range of 
magnitudes with a median I = 22.1. The present deeper I- 
band observations remain consistent with this result. The 
median I for the HDF-N and SSA13 1.4 GHz sources with 
/-band measurements is 22.37'fQ'\\ (68% conhdence). The 
distribution implies that the bulk of the radio sources (just 
over 60%) should be bright enough (/ ^ 23.5) for spectro¬ 
scopic identihcation. The remaining sources will mostly be 
too optically faint for spectroscopic identihcation and will 
require another method for determining their redshifts. 
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Fig. 2.— Redshift versus 1.4 GHz flux for the spectroscopically- 
observed radio sources in the HDF-N and SSA13 flelds. The solid 
squares denote sources with redshift identifications, and the solid 
triangles at 2 : = 7 denote sources for which no redshifts could be 
obtained. The large open squares denote the median redshifts for 
the identified sources. The 68% confidence range on the median red¬ 
shifts are indicated by the error bars. The large solid square with 
error bars shows the median redshift and 68% confidence range of 
the millijansky radio sample of Waddington et al. (2001). 

In Fig. 2 we plot redshift versus S'i, 4 ghz for the 
spectroscopically-observed radio sample. Sources with 
redshifts are denoted by solid squares, and sources for 
which no redshifts could be obtained are denoted by solid 
triangles at z = 7. The large open squares denote the 
median redshifts for the identified sources. The median 
redshift distribution is essentially flat across the entire mi¬ 
cro jansky regime probed. We also plot the median red¬ 
shift (large solid square) from a millijansky radio study 
by Waddington et al. (2001). This point shows a con¬ 
tinuation of the flat redshift distribution to brighter radio 
fluxes. The selection of optically-bright sources inherent in 
a spectroscopically-identified sample means we are likely 
missing the high redshift tail of the distribution function. 
This tail could be skewing to higher redshifts as we move 
to fainter radio fluxes while the main body of the redshift 
distribution remains invariant. We shall return to this 
point when we discuss the submillimeter properties of the 
radio sources. 

In Fig. 3 we plot rest-frame absolute /-band mag¬ 
nitude, M/, calculated using a standard spiral galaxy 
//-correction, versus redshift for the spectroscopically- 
identifled radio sources (solid squares). The effect of the 
/ ~ 23.5 spectroscopic limit on Mj is indicated by the 
dashed line, and the effect of our growing incompleteness 
beyond / ~ 22 is indicated by the dotted line. Fig. 3 
shows a remarkably narrow range in M/ with redshift (as 
indicated by the solid lines) that suggests the radio sources 
are chosen from approximately L* optical galaxies, regard¬ 
less of redshift, and the radio hosts are, at a crude level, 
standard candles in their optical properties. In the red¬ 
shift range z = 0.1-0.5 the mean Mj = —22.0 ± 0.9; for 
2 = 0.5 — 0.8 the mean Mj = —22.6 ± 0.8 (by contrast, 
classical radio galaxies are approximately 5-10 xL* optical 
galaxies). As can be seen from Fig. 3, the small increase 
in absolute magnitude at higher redshifts largely may be 
due to the apparent magnitude selection effects. 

To verify that this result is not an effect of our flux 


limited radio sample, we construct a toy model whereby a 
Schechter luminosity function with Mj = —22.5 is evolved 
in luminosity, increasing with redshift as (l-|-z)^ from 
2 = 0.1 — 1.5, and normalized by the volume element, dV. 
If we were to assume that radio luminosity is proportional 
to optical (/-band) luminosity, our observed range in Mj 
would result in only sub-T* radio sources lying in our sam¬ 
ple. However, if we apply our radio flux limit (~40y^Jy) to 
the model, we notice that we can detect only those radio 
sources with luminosities greater than L* at all redshifts 
> 0.3, with 70% of the detectable sources in the observed 
volume lying outside of the equivalent Mj range. We con¬ 
clude that our observed result (micro) ansky radio sources 
being hosted by similarly luminous optical galaxies) is not 
a selection effect of our flux limited radio sample. 



Redshift 

Fig. 3.— Mi versus redshift for the radio sources in the HDF-N 
and SSA13 fields with spectroscopic redshifts (solid squares). The 
solid lines show the approximate range in Mj from —21 to —24 
where the majority of the sources lie. The dashed line illustrates 
the effect of an 7 ~ 23.5 spectroscopic limit on Mj, and the dotted 
line illustrates that of an 7 22 selection. 

It is not clear why the radio sources should all lie in high 
optical luminosity galaxies. At the higher redshifts the re¬ 
stricted absolute magnitude range can be understood from 
the selection effects: the lower cut-off is due to low optical 
luminosity galaxies being faint and not spectroscopically 
identified, and the upper cut-off reflects the small number 
of galaxies with optical magnitudes above the optical L*. 
However, at the lower redshifts there is no selection against 
low optical luminosity galaxies, and the absence of these 
galaxies is striking. Many of the low optical luminosity 
galaxies at low redshifts have blue colors and substantial 
ongoing star formation that might be expected to produce 
substantial synchrotron emission; thus, we might expect 
this population to correspond to a large part of the ra¬ 
dio population. The absence of these objects argues that 
they are weak in the radio relative to their optical (and 
ultraviolet) luminosities. 

This presumably is the same effect that is seen in the 
local far-infrared (FIR)-radio correlation, which is not pre¬ 
cisely linear; low optical luminosity galaxies have lower 
radio luminosities than expected (Fitt et al. 1988; Cox et 
al. 1988; Devereux & Bales 1989). One likely explanation 
is that small galaxies are not very effective radio sources 
because cosmic rays are more likely to escape by diffu¬ 
sion or convection (Chi & Wolfendale 1990; Condon et al. 
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1991.) As our sources are selected from a flux-limited radio 
survey, low optical luminosity galaxies will be eliminated 
preferentially from our sample. This is a serious problem 
for radio estimates of the local star formation rate den¬ 
sity (e.g., Mobasher et al. 1999; Haarsma et al. 2000), as 
a substantial fraction of the ultraviolet luminosity density 
is generated by faint sub-L* galaxies at low redshifts (e.g., 
Cowie, Songaila, & Barger 1999; Wilson et al. 2002). The 
radio-based star formation rate density estimates already 
tend to lie at the high end of the values measured in the 
z = 0-1 range, but despite this, the lowest redshift values 
must be substantially underestimated. 

4. SUBMILLIMETER PROPERTIES OF THE MICROJANSKY 
RADIO SOURCES 

In Fig. 4 we plot Sgso/im versus J-band magnitude for 
the radio sources with submillimeter measurements, and 
in Fig. 5 we plot S'sso/im versus 51.4 ghz. The data are 
divided both according to observation technique (trian¬ 
gles for jiggle map or raster mode, circles for photometry 
mode) and according to submillimeter detection signifi¬ 
cance (large solid symbols for > 3cr, small open symbols 
for < Str). Analyses of submillimeter source reliability 
from SCUBA maps suggest a threshold of 4cr (Bales et al. 
2000; Scott et al. 2002; Borys et al. 2002; Cowie, Barger, 
& Kneib 2002); thus, we have enclose our > 4cr sources 
in a second larger symbol. However, confidence in a tar¬ 
geted submillimeter detection of a known optically-faint 
radio source is much higher than in a randomly detected 
submillimeter source. For the 278 radio sources with sub¬ 
millimeter measurements, we expect only 0.5 spurious 3cr 
sources. 



Fig. 4. — Submillimeter flux versus /-band magnitude for the 
microjansky radio sources in the HDF-N, SSA13, and SSA22 fields 
with submillimeter measurements. Circles denote SCUBA photom¬ 
etry observations and triangles denote SCUBA jiggle map or raster 
map observations. Large filled symbols denote sources with > Sc 
submillimeter fluxes. Sources with > 4 (t submillimeter fluxes are 
enclosed in a second larger symbol. The error-weighted mean sub¬ 
millimeter fluxes are shown as large open squares with Icr uncer¬ 
tainties. The vertical dashed line at / = 23.5 shows our division be¬ 
tween optically-faint and optically-bright sources. Error bars have 
not been plotted but are typically ~ 2 mjy. 


4.1. Statistical analyses 

In Fig. 4 we show the error-weighted mean submillime¬ 
ter fluxes (open squares with Icr uncertainties) for five 
bins of 53 sources at successively fainter /-band mag¬ 
nitudes. The sample as a whole has an error-weighted 
mean 5850/jm= 1-73 ± 0.09 mJy. Since the distribution of 
the noise is not Gaussian (see Cowie, Barger, & Kneib 
2002 ), the computed error may underestimate the true 
value. However, when Monte Carlo simulations of an 
equal number of randomly distributed sources are ana¬ 
lyzed in the same way, the high significance of the re¬ 
sult is confirmed. The optically-faint (/ > 23.5) sub¬ 
sample is highly significant at 2.37 ± 0.13 mJy, while the 
optically-bright (/ < 23.5) subsample is also significant at 
1.09±0.13mJy. Thus, the optically-faint radio sources are 
about two times more submillimeter luminous, on average, 
than the optically-bright radio sources. While we have 
somewhat oversampled the optically-faint radio sources 
with our targeted submillimeter photometry observations, 
the result is consistent when we consider only our contigu¬ 
ous submillimeter jiggle map regions. 

Both the optically-faint and optically-bright radio- 
selected samples have higher submillimeter fluxes than X- 
ray-selected samples: for the 1 Ms Chandra Deep Field- 
North X-ray sample the error-weighted mean 5850/im= 
1.21 ± 0.14mJy (Barger et al. 2001b), and for the 75 ks 
Elais-N2 sample, 1.25±0.40mJy (Almaini et al. 2002). In 
fact, Barger et al. (2001b) found that the submillimeter 
signal of their X-ray sample was largely driven by the sub¬ 
set of sources that overlapped with optically-faint radio 
sources. 



®1.4 GHz (/^Jy) 


Fig. 5. — Submillimeter flux versus radio flux for the microjansky 
radio sources in the HDF-N, SSA13, and SSA22 fields with submil¬ 
limeter measurements. Symbols are as in Fig. 4. 

In Fig. 5 we show the error-weighted mean submillime¬ 
ter fluxes (open squares with la uncertainties) at suc¬ 
cessively fainter 1.4 GHz fluxes. The mean rises (from 
5850 Aim~ 1 mJy to 3mJy) with increasing radio flux out 
to about 100/iJy before dropping to 5850/rm=l-2 mJy for 
the brightest radio sources. For radio fluxes approaching 
the millijansky regime, we expect an increasing contribu- 
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tion from active galactic nuclei (AGN). This appears to 
be reflected in the trend of the means: the brightest radio 
sources in our sample are, on average, not the strongest 
submillimeter emitters (luminous AGN typically have low 
5'850/im/5'i,4GHz ratios). Since the SSA13 radio image is 
deeper than the HDF-N and SSA22 images, the number 
of radio sources probed in the submillimeter is not sim¬ 
ply proportional to the radio source counts. However, the 
statistics should not be affected by this, unless significant 
large-scale structure dominates the small field covered to 
this depth. Note that within the area probed by our sub¬ 
millimeter observations, only two radio sources have radio 
fluxes higher than 1 mJy. 

4.2. Direct submillimeter source detections 

Fifty of the radio sources in our sample are significantly 
(> 3cr) detected individually in the submillimeter. Of 
these, 38 fall into the optically-faint (/ > 23.5) region of 
Fig. 4, which is consistent with the conclusion of previous 
SGUBA surveys that most submillimeter sources are opti¬ 
cally faint (e.g., BGROO; Smail et al. 2002; GOl). The per¬ 
centage of 5850/im> 5mJy (> 3 (t) submillimeter sources 
detected in our survey versus optical brightness is pre¬ 
sented in Table 1; the radio selection technique becomes 
highly efficient for optical magnitudes fainter than I ~ 24. 

While the submillimeter detection rate of / < 23.5 ra¬ 
dio sources is rather poor, this population still makes up a 
fraction of our submillimeter sample. To quantify this, we 
consider only regions of our survey that were mapped with 
SGUBA to avoid the bias introduced by the optically-faint 
criterion used for the targeted photometry follow-up. Us¬ 
ing our combined jiggle and raster maps of the three fields, 
we find that 20 of the submillimeter-detected sources have 
I > 23.5 and 9 have I < 23.5. Thus, 69±9%, where the er¬ 
rors correspond to the 1 sigma variance in the expectation 
value, of the radio-submillimeter population has / > 23.5. 
This appears to be inconsistent with the claim by Ivison et 
al. (2002), based on a smaller sample, that the bulk of the 
radio-detected submillimeter sources are optically bright. 

4.3. Overlap of the bright submillimeter population with 
the microjansky radio population 

Using our jiggle and raster mapped regions (> 
300 arcmin^), we can determine what fraction of bright 
submillimeter sources are picked out by an ultradeep ra¬ 
dio sample. As the radio properties in the SSA22 field are 
highly non-uniform due to bright radio sources, we again 
restrict our analysis to the HDF-N and SSA13 fields. We 
also restrict to > 4tT sources with Ss 50 tJ.m> 5mJy. We 
have adopted an 8" search radius to avoid missing asso¬ 
ciations due to the positional uncertainties of the S/N^4 
SGUBA detections. Note that the 8" search radius should 
not give rise to any spurious associations as the surface 
density of radio sources is only ~ 2 arcmin“^ in the cen¬ 
tral regions of our maps. We find that 66 ± 7% of bright 
submillimeter sources (29 of the /N~4 SGUBA detections) 
have radio identifications. This is consistent with Ivison 
et al. (2002), who found that 18 of the 30 submillimeter 
sources in their sample have radio identifications, which 
translates to a percentage of 60 ± 9%. From § 4.2, about 
70% of these counts will be picked out by the optically- 
faint radio sources. 



Fig. 6.— I — K versus the radio of submillimeter to radio flux 
for the microjansky radio sources in the HDF-N and SSA22 fields 
with both submillimeter and /C-band measurements. Sources with 
marginal /C-band detections are shown as 3cr upper limits. The 
dashed line shows the Very Red Object limit, I — K = 3.5. 

4.4. I — K Colors 

We have AT-band measurements for a subset of our radio 
sample. Both the optically-bright and optically-faint radio 
subsamples have significantly redder I — K colors than the 
field galaxy population. The optically-faint radio sources 
show a monotonic trend to redder colors (Richards et al. 
1999; BGROO): at / = 22, median I — K = 2.9 and at 
/ = 24.5, median I — K = 3.4. For comparison, the field 
galaxy sample has median / —A'~2 at/~22 and median 
/ — AT ^ 2.3 at / ~ 24.5 (e.g., Barger et al. 1999). 

In Fig. 6 we plot I—K versus 5850/im/5i.4GHz for the 3a 
submillimeter-detected radio sources with A'-band mea¬ 
surements. All of the sources have J-band detections, but 
three of the sources are only marginally detected at K and 
are shown as 3cr upper limits. The submillimeter-detected 
sources have median I — K = 4.3 and median K = 20.7. 
From Fig. 6 we can see that a color cut of / — AT >3.5 
includes 95% of the submillimeter-detected sources. This 
is consistent with the results of Wehner, Barger, & Kneib 
(2002), who used a submillimeter survey of three massive 
cluster fields to determine that the median submillimeter 
source has K = 20 — 20.3, depending on the cluster sample 
used, and that I—K > 3.5 sources (Very Red Objects) con¬ 
tribute about two-thirds of the 850 ^mextragalactic back¬ 
ground light. 

There is no obvious trend of / — AT color with 
^850 pm /Si. 4 GHz in Fig. 6, implying that the near-infrared 
colors of the submillimeter-selected galaxies are insensitive 
to redshift over the z = 1-3 range as would be expected if 
the red colors are primarily caused by dust extinction. 

5. MILLIMETRIC REDSHIFT ESTIMATION FOR THE 
MICROJANSKY RADIO-SELECTED SAMPLE 

5.1. Evolution of the radio population with redshift 

The radio power of the spectroscopically-identified sam¬ 
ple in the HDF-N and SSA13 fields is shown in Fig. 7, 
where we have computed the radio AT-correction assuming 
a synchrotron spectral index of —0.8. The selection limits 
of the HDF-N and SSA13 samples (40 ^Jy and 25 /rJy, re¬ 
spectively) are shown by the solid curves. The conversion 
to FIR luminosity, computed using the local FIR-radio cor¬ 
relation (e.g., Gondon 1992), is shown on the right-hand 
axis, and the dotted and dashed lines show the FIR lu¬ 
minosity limits for luminous infrared galaxies (LIGs) and 
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ultraluminous infrared galaxies (ULIGs), respectively. 

There is a rapid evolution in the upper bound of the ra¬ 
dio power. This is the well-known effect that large galaxies 
at z ~ 1 frequently have very large star formation rates 
while at lower redshifts the rates turn down rapidly (Cowie 
et al. 1996, where the phenomenon is called downsizing). 
An equivalent description is that near-ULIGs or ULIGs 
are extremely rare at low redshifts, but the number den¬ 
sity rises rapidly with redshift so that they are relatively 
common by z = 1 (e.g., Blain et al. 1999; BGROO). In 
combination with the lower-bound selection imposed by 
the flux limits, this effect means that we have a rapidly 
changing but rather well-defined population with redshift. 
At low redshifts (z < 0.4) we are selecting Milky Way type 
galaxies, at intermediate redshifts we are selecting LIGs, 
and at z > 1 we are selecting only near-ULIGs or ULIGs. 



Fig. 7. — Radio power versus redshift for the radio sources in the 
HDF-N and SSA13 fields with spectroscopic redshifts. The solid 
curves show the selection limits for the two fields (40 /iJy and 25 /J-Jy, 
respectively). The right vertical axis shows the logarithm of the FIR 
luminosity (in solar luminosities) that would be inferred from the 
radio power, based on the local FIR-radio correlation. The FIR lu¬ 
minosity limits of LIGs (10^^ Lq) and ULIGs (10^^ Lq) are shown 
by the dotted and dashed lines, respectively. 

5.2. Average Ssbo fi,m/S iagrz flux ratios 

Our present dataset enables us to determine the nor¬ 
malization with redshift of the error-weighted mean 
'S' 850 /im/'S'i. 4 GHz Aux ratio. In Fig. 8 we plot the loga¬ 
rithm of the error-weighted mean S' 850 /im/<S'i. 4 GHz for the 
radio sources with spectroscopic redshifts and submillime¬ 
ter measurements (large open squares). The four > Scr 
submillimeter-detected sources with redshifts (see § 5.3) 
are shown as dots. We include on the figure the tracks with 
redshift for three spectral energy distributions (SEDs) 
with different dust properties: Arp 220 (short dashed), 
a typical lUE starburst (solid), and a Milky Way type 
galaxy (dotted). We have calculated these tracks using the 
Dale et al. (2001) and Dale & Helou (2002) empirical SED 
library. We also include three tracks from the literature 
(DCEOO; GYOO; Yun & Carilli 2002) that were derived us¬ 
ing different galaxy samples and assumptions; these tracks 
are the most widely used in millimetric redshift estimation. 
As would be expected from the luminosity selection, the 
error-weighted means follow the track of a quiescent to 
moderately star-forming galaxy at lower redshifts, while 
at higher redshifts the observations migrate to values Icr 


below the Arp 220 and GYOO tracks. (We note that Yun, 
Reddy, & Gondon (2001) found that sources like Mrk 231, 
whose AGN has a radio luminosity comparable to its star- 
burst luminosity, likely constitute a few to 10% of the sub¬ 
millimeter population. Such sources may lie preferentially 
in our radio-selected sample and hence have lowered our 
high-redshift ULIG points by about 0.3 in dex.) 



redshift 

Fig. 8. — log 58504GHz versus redshift for the microjansky 
radio sources with spectroscopic redshifts and submillimeter mea¬ 
surements. Squares denote the error-weighted means in six equal- 
number bins. The four > Sc submillimeter-detected sources with 
reliable redshifts are shown as dots. The tracks are Arp 220, a Milky 
Way type galaxy, a typical lUE starburst, and Dunne et al. (2000), 
Carilli & Yun (2000), and Yun &; Carilli (2002). 

The radio power and the FIR luminosity correlate lo¬ 
cally with dust temperature, and there is no evidence that 
this relation varies significantly with redshift (Chapman 
et al. 2003, Blain, Barnard & Chapman 2003). Given 
this, the evolution in radio power and in S'850/im/5'i,4GHz 
with redshift show that for the microjansky radio popu¬ 
lation, millimetric redshift estimates at low redshifts are 
best made with a template intermediate between a Milky 
Way type template and a starburst template, and those at 
high redshifts (z > 1) with an Arp 220 template. Thus, 
we can apply the millimetric redshift technique to other 
radio-submillimeter populations provided we have knowl¬ 
edge of the characteristic FIR luminosity or redshift of the 
population before we choose a template SED. 

In Fig. 3 we showed that nearly all radio sources with 
redshifts have similar M/. This suggests that the optically- 
faint radio sources, rather than being low optical luminos¬ 
ity galaxies at moderate redshifts are instead luminous 
galaxies at z > 1, where the negative AT-correction in 
the submillimeter provides a nearly one-to-one mapping 
of submillimeter flux to submillimeter luminosity. For the 
sources that are individually detected in the submillimeter, 
the 5850^1111 fluxes (see Fig. 4) give FIR luminosities that 
are at least as large as Arp 220 (> 10^^ Lq). We therefore 
can use the Arp 220 template in making millimetric red¬ 
shift estimates for the optically-faint microjansky sources. 

In Fig. 9a we show S' 850 /im/<S'i. 4 GHz versus I.4GHz flux. 
The right-hand vertical axis shows the redshift that would 
be inferred from the Arp 220 template. In Fig. 9b we 




















divide the sample into optically-bright and optically-faint 
subsamples. We can see from this figure that the optically- 
faint subsample does indeed move to higher redshift as we 
probe to fainter radio fluxes. 

100 1000 





Fig. 9. — (a) Ss50fim/SiAGHz versus 1.4GHz flux for the micro- 
jansky radio sources in the HDF-N, SSA13, and SSA22 fields with 
submillimeter measurements. Large filled symbols denote sources 
with > 3(J submillimeter fluxes. Triangles (circles) denote sources 
with I < 23.5 (7 > 23.5). Assuming the sources span a rea¬ 
sonably small range in dust temperature and properties, the y- 
axis can be interpreted as a redshift measure (right-hand scale). 
The redshift range covered by our radio sources is 2 0 — 3 for 

typical dust temperatures encountered in local starburst galaxies 
and ULIGs. The mean 5'850/2in/5'i.4GHz is independent of radio 
flux for sources fainter than about 200 pJy, within Icr statistical 
errors, (b) Mean 5850Aim/5'i.4 ghz versus 1.4GHz flux for the 
optically-bright and optically-faint subsamples. The nearly con¬ 
stant 5850 ^m/5i.4 GHz ratio of the total sample is dominated by 
the optically-bright sources. 


Chandra Deep Field-North exposure (Barger et al. 2001b). 
The spectroscopic redshifts of this and the z = 0.510 
source are lower than their millimetric redshift estimates 
(assuming an Arp 220 SED, we find z = 1.6 and z = 1.9, 
respectively). Sources with dust temperatures cooler than 
the Arp 220 ~ 45 K template will have true redshifts and 
FIR luminosities lower than predicted (see also Yun & Car- 
illi 2002). The z = 1.013 galaxy exhibits a cool dust tem¬ 
perature (~30 K), similar to recently identified cold ULIGs 
from the FIRBACK 170/rm ISOPHOT survey (Chapman 
et al. 2002c). In the case of the z = 0.510 galaxy, however, 
the dust would have to be prohibitively cold (~ 18 K) for 
the millimetric redshift to match the spectroscopic identi¬ 
fication. This source may therefore be a misidentification, 
with the true source being lensed by the bright foreground 
galaxy (Chapman et al. 2002b). 


0.1 0.5 1 



redshift 


Fig. 10.— Submillimeter flux versus redshift for the microjan- 
sky radio sources with submillimeter measurements. The I < 23.5 
sources with (without) spectroscopic redshifts are denoted by solid 
(open) triangles, and the I > 23.5 sources with submillimeter detec¬ 
tions (and hence millimetric redshifts) are denoted by solid circles. 
The sources without > Sa submillimeter detections are denoted by 
left-pointing arrows at the millimetric redshift limit (estimated us¬ 
ing the 2cr submillimeter uncertainty) and the measured 5850 
flux. The error-weighted mean 5850Aiin/5i.4GHz of these sources is 
denoted by a large open circle at its millimetric redshift estimate of 
2 = 1.13. The track of a galaxy twice as luminous as Arp 220 is 
shown as a dashed line. The expected characteristic error at 2 : = 2, 
based on the range in properties for local ULIGs relative to the 
Arp 220 template, is indicated by a horizontal bar. 


5 . 3 . Individual sources 

There are four significant (> 3cr) submillimeter sources 
in our sample with spectroscopic redshifts (z = 3.405, 
z = 2.565, z = 1.013, and z = 0.510); all have I < 23.5. 
The z = 3.405 source has Lyman a and CIV in emission, 
and the z = 2.565 source is a broad-line quasar that was 
detected in the 100 ks Chandra exposure of the SSA13 
field (Barger et al. 2001a). The redshifts of both sources 
are very close to their millimetric redshift estimates. The 
z = 1.013 source shows strong Balmer absorption features 
and [Ne III] in emission and was detected in the 1 Ms 


5.4. Combined spectroscopic and millimetric redshift 

sample 

In Fig. 10 we plot Sssofim versus spectroscopic (solid tri¬ 
angles) or millimetric (solid circles) redshift for the sources 
with submillimeter measurements. We also show a rough 
error bar based on the range in properties for local ULIGs 
relative to the Arp 220 template. We have very little infor¬ 
mation on the spectroscopically-unidentified radio sources 
that are not individually detected in the submillimeter. 
We can put conservative limits on their redshifts by esti- 
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mating millimetric redshifts from the CYOO relation and 
the 2a submillimeter uncertainty. We denote these sources 
by the left-pointing arrows in Fig. 10. 

We can go one step further by determining a single av¬ 
erage point for this marginally detected population from 
the ratio of the error-weighted mean submillimeter flux 
and the error-weighted mean radio flux. We denote this 
point by a large open circle a,t z = 1.13 on Fig. 10. As we 
have already removed all the Ss 50 ii.m> 3ct sources, submil¬ 
limeter sources without radio counterparts are not present. 
Since the steep submillimeter source counts imply a large 
Eddington bias, it is important not to overinterpret this 
result, but we note that it falls naturally within the dis¬ 
tribution of the more robustly identified sources. 

Comparing the data with the track followed by a ULIG 
with twice the luminosity of Arp 220 (dashed line in 
Fig. 10), we can see the rapid rise in luminous FIR sources 
with increasing redshift. Below z ~ 1 there are essentially 
no such objects, while above z ~ 1 they are common. 

SUMMARY 

We presented optical and submillimeter (jiggle and 
raster maps and photometry) imaging and optical spec¬ 
troscopy for a sample of microjansky radio sources de¬ 
tected in ultradeep 1.4 GHz maps of the HDF-N, SSA13, 
and SSA22 fields. Our spectroscopic identifications for 169 
of the sources in the HDF-N and SSA13 fields revealed 
a flat median redshift distribution across the entire mi¬ 
crojansky regime probed. The identified sources show a 
remarkably narrow range in rest-frame absolute J-band 
magnitude that suggests they are chosen from approxi¬ 
mately L* optical galaxies, regardless of redshift. At the 
higher redshifts the narrow range can be understood from 
selection effects, but at the lower redshifts there is no se¬ 
lection against low optical luminosity galaxies. Thus, the 
absence of low optical luminosity galaxies in the sample ar¬ 
gues that they are weak in the radio relative to their optical 
and ultraviolet luminosities, possibly due to small galaxies 
not being efficient at retaining the cosmic rays necessary 
to host microjansky radio sources. The fact that these 
galaxies are eliminated preferentially from a flux-limited 
radio survey poses a serious problem for radio estimates 
of the local star formation rate density since a substantial 
fraction of the ultraviolet luminosity density is generated 
by sub-L* galaxies at low redshifts. 

From a statistical analysis of our submillimeter measure¬ 
ments for 278 of the microjansky radio sources in the HDF- 
N, SSA13, and SSA22 fields, we found that optically-faint 
radio sources are about three times more submillimeter 
luminous, on average, than optically-bright radio sources, 
but both populations are, on average, very significantly 
detected in the submillimeter. Fifty of the radio sources 
in our sample are significantly (> 3 ct) detected individu¬ 
ally in the submillimeter; 38 of these are optically faint 
(/ >23.5). We have spectroscopic redshifts for four of the 
I < 23.5 submillimeter-detected sources, three of which 
are in the range z = 1-3.4; the fourth source is at z = 0.510 
and may be a misidentification. 

To quantify the percentage of optically-bright and 
optically-faint radio sources that are significantly detected 
in the submillimeter, we considered only the regions of 
our survey that were mapped with SCUBA (as opposed 
to the regions that were targeted with photometry). We 


found that 69 ± 9% of the radio-submillimeter popula¬ 
tion has I > 23.5. We also found that 66 ± 7% of the 
Ss 50 fj,ni> 5mJy (> dcr) sources in our fields are radio- 
identified. 

We used our spectroscopically-identified microjansky ra¬ 
dio sample to determine the evolution with redshift in the 
radio power and hence the FIR luminosity (through the 
FIR-radio correlation). We found that ULIGs or near- 
ULIGs are extremely rare at low redshifts but are rela¬ 
tively common by z = 1. In combination with the lower- 
bound selection effects, this means that we have a pop¬ 
ulation which is rapidly changing with redshift but that 
is fairly tightly selected at each redshift. At low redshifts 
(z < 0.4) we are selecting Milky Way type galaxies, at in¬ 
termediate redshifts we are selecting LIGs, and at high red¬ 
shifts (z > 1) we are selecting only ULIGs or near-ULIGs. 
From the error-weighted mean 5'850/jm/>5'i.4GHz hux ra¬ 
tio we see the same evolution with redshift, such that at 
lower redshifts the error-weighted means follow the track 
of a quiescent to moderately star-forming galaxy, while at 
higher redshifts the observations migrate to values consis¬ 
tent with Arp 220. Thus, millimetric redshift estimates 
at low redshifts are best made with a template intermedi¬ 
ate between a Milky Way type template and a starburst 
template, while at high redshifts (z > 1) an Arp 220 tem¬ 
plate is appropriate. With this knowledge, we estimated 
millimetric redshifts or limits for all of the optically-faint 
microjansky radio sources in our sample with submillime¬ 
ter measurements. 
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Table 1 

Average Submillimeter Fluxes for the HDF-N+SSA13+SSA22 Sample in /-magnitude bins 


7-mag Bin 

5*850 /tm (mJy) 

CT850 (mjy) 

number in bin 

% SMM detections (> 30-) 

7 > 25 

2.38 

0.18 

68 

26% 

24 < 7 < 25 

2.85 

0.23 

43 

37% 

23 < 7 < 24 

1.32 

0.27 

32 

24% 

22 < 7 < 23 

1.00 

0.28 

36 

8% 

21 < 7 < 22 

1.13 

0.26 

31 

9% 

19.5 < 7 < 21 

1.10 

0.31 

31 

3% 

7 < 19.5 

1.08 

0.29 

37 

5% 



